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Abstract. Pulsar timing experiments are reaching sufficient sensitivity to detect 
a postulated stochastic gravitational wave background generated by merging 
r- ' supermassive black hole systems in the cores of galaxies. We describe the techniques 

<*S |. behind the pulsar timing detection method, provide current upper bounds on the 

q \ amplitude of any gravitational wave background, describe theoretical models predicting 

%-{ ' the existence of such a background and highlight new techniques for providing a 

statistically rigorous detection of the background. 



> 

' 1. Introduction 

O ' 

It is now possible to make timing observations of millisecond pulsars to a precision of 
! ~ 100 ns. Such exquisite timing precision allows the pulsar astrometric, spin and orbital 
parameters to be determined with great accuracy. This has recently allowed us to limit 
O ! the rate of change of the gravitational constant, \G/G\ < 2.1 x 10~ 12 yr™ 1 , and to obtain 
the most accurate astrophysical distance estimate outside of the Solar System (Verbiest 
et al. 2007, submitted to ApJ). Multi- frequency pulsar observations have allowed us 
to study the Solar corona (You et al. 2007, in press) and density fluctuations in the 
interstellar plasma (You et al. 2007b). One of the most exciting applications of such 
data-sets is to search for the signatures of gravitational waves (GWs) passing over the 
Earth. This is the main goal of the Parkes Pulsar Timing Array (PPTA) project^, which 
aims to observe 20 millisecond pulsars with a timing precision close to 100 ns over more 
than five years. 

The existence of a GW signal can be identified by analysis of pulsar timing residuals. 
In brief, the pulsar timing model (see, e.g. Edwards, Hobbs & Manchester 2006) allows 
the observed arrival times of pulses from a radio pulsar to be compared with a model of 
the pulsar's astrometric, orbital and rotational parameters. The differences between the 
actual arrival times and the predicted arrival times are the "pulsar timing residuals". 
Any features observed in the timing residuals indicate the presence of unmodelled 
effects which may include calibration errors, (additional) orbital companions, spin-down 



J http://www.atnf.csiro.au/research/pulsar/ppta 
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Figure 1. Characteristic strain sensitivity for existing and proposed GW detectors as 
a function of GW frequency. Predicted signal levels from various astrophysical sources 
are shown. 

irregularities, or GWs. The induced pulsar timing residuals due to a GW signal were 
first calculated by Sazhin (1978) and Detweiler (1979). They showed that a GW signal 
causes a fluctuation in the observed pulse frequency bvjv which affects pulsar timing 
residuals at time t from the initial observation as 

m = - j^dt. a) 

For a given pulsar and GW source, this effect is only dependent upon the characteristic 
strain at the pulsar and at the Earth. The GW strains evaluated at the positions of 
multiple pulsars will be uncorrelated, whereas the component at the Earth will lead to a 
correlated signal in the timing residuals of all pulsars. As pulsars are typically observed 
every few weeks for many years and standard pulsar timing techniques absorb any 
low-frequency GWs by fitting for the pulsar's spin-dowrj§], pulsar timing experiments 
are sensitive to GW signals in the ultra low frequency (/ ~ 1CT 9 Hz) band. This 
technique is therefore complementary to other GW detection methods such as the Laser 
Interferometer Space Antenna (LISA) and ground based interferometer systems (such 
as LIGO and VIRGO), which are sensitive to higher frequency GWs (/ ~ 10~ 2 ' 5 and 
f ~ 10 2 5 Hz respectively). In Figure [1] we plot the sensitivity of the pulsar timing 
experiments, LISA and LIGO and the expected GW sources. For the pulsar timing 
experiments, maximum sensitivity is reached at / = 1/T where T is the data span. As 
described above, the pulsar experiments are not sensitive to GWs with f < 1/T and, 
for white pulsar timing residuals, the upper bound on the characteristic strain increases 

§ As the intrinsic pulsar spin-period and rate of spin-down is unknown, a quadratic curve is always 
fitted and subsequently removed from the pulsar timing residuals. 
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Figure 2. Simulation of the induced timing residuals for PSR B1855+09 due to the 
postulated binary black hole system in the radio galaxy 3C66B. 

as f 3/2 for / > 1/T. 

In this paper we review the potential sources of GWs (§2]) before describing how 
upper bounds on the amplitude of any existing background can be calculated and 
highlighting the astrophysical implications of these bounds. In £J3] we describe how a 
stochastic GW background may be detected. Finally, we highlight some practical issues 
that need to be addressed in achieving our sensitivity goal and describe possible 
future experiments ($5]). 

2. Potential sources of gravitational waves 

Theoretical models predict that pulsar timing experiments will be sensitivity to burst 
GW sources, individual sources and to a stochastic background. 

2.1. Single sources 

Recently Sudou et al. (2003) reported the possible discovery of a binary supermassive 
black hole system in the radio galaxy 3C66B. As shown by Jenet et al. (2004), the 
GW radiation from this postulated system would have produced clearly detectable 
fluctuations in the timing residuals of existing pulsar data-sets. The non-detection of the 
expected GW signal (see Figure [2]) ruled out the postulated system with 95% confidence. 
However, detecting (or limiting the existence of) individual supermassive binary black 
hole systems is necessary for testing theories of hierarchical galaxy formation. Order- 
of-magnitude calculations demonstrate that a system with a chirp mass of 10 9 M Q and 
an orbital period of 10 yr, at a distance of 20Mpc will induce a sinusoidal signal in the 
timing residuals with amplitude ~ 100 ns at a frequency of twice the orbital frequency. 

2.2. Burst sources 

Sources of burst GW emission that may be detectable with current and proposed 
pulsar timing experiments include 1) the formation of supermassive black holes (Thorne 
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Figure 3. The solid line indicates our current limits on GW backgrounds of specified 
spectral exponent. The dashed line underneath indicates the expected limits within 
5 years. The near vertical lines give predicted background amplitudes. The "star"- 
symbol indicates the earlier limit obtained by Kaspi et al. (1994). 



& Braginskii 1976), which lead to a day-long burst of radiation, 2) highly eccentric 
supermassive black hole binaries (Enoki & Nagashima 2007), 3) close encounters of 
massive objects (Kocsis et al. 2006) and 4) cosmic string cusps (Damour & Vilenkin 
2001). The sensitivity of a pulsar timing array experiment to burst GW emission 
depends upon the number of pulsars in the array, the timing precision, number of 
observations, the pulsar positions on the sky and the coordinates of the GW burst source. 
With our final data-sets we expect to achieve an angular resolution of approximately 30° 
(Lommen, private communication). By sharing data with timing array projects being 
carried out in the Northern hemisphere we expect to improve both the angular resolution 
of our detector and our sensitivity to burst sources in the Northern Hemisphere. 



2.3. Stochastic background 

A stochastic background of GWs due to binary supermassive black holes (Jaffe & Backer 
2003, Wyithe & Loeb 2003, Enoki et al. 2004, Sesana et al. 2004), cosmic strings or 
relic GWs from the big bang (Maggiore 2000) are all potentially detectable. In most 
models, the GW strain spectrum, h c (f), can be represented by a power-law in the GW 
frequency, /, 

where the spectral exponent a = —2/3, —1 and —7/6 for likely GW backgrounds due 
to coalescing black hole binaries, cosmic strings and relic GWs respectively (Figure [3]). 
The energy density of the background per unit logarithmic frequency interval can be 
written as 

*W/) = |^Ac(/) 2 (3) 

3 H 

where H is the Hubble constant. 
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Jenet et al. (2006) introduced a method to place an upper bound on A for a 
given a using actual pulsar timing residuals. In contrast with earlier techniques (e.g. 
Kaspi et al. 1994) this method used the tempo2 software package (Hobbs, Edwards 
& Manchester 2006) to account for all the fitting procedures undertaken during the 
pulsar timing process. Using observations from the PPTA project along with archival 
data from the Arecibo telescope, Jenet et al. (2006) obtained an upper limit on the 
energy density per unit logarithmic frequency interval of flQ^[l/(8jT)]h 2 < 2 x 10 -8 for 
a = —2/3, which corresponds to A < 10~ 14 (see Figure [3] and the "Current limit" in 
Figured]). This is the most stringent limit to date on the existence of a GW background 
in the nano-Hertz frequency band and has astrophysical and cosmological implications. 
For instance, this limit was used to 1) constrain the merger rate of massive black hole 
systems at high redshift, 2) rule out some relationships between the black hole mass 
and galactic halo mass, 3) constrain the rate of expansion in the inflationary era and 4) 
provide an upper bound on the dimensionless tension of cosmic strings. 

We have recently shown that the Jenet et al. (2006) technique described above 
suffers from two limitations. First, it does not give the lowest possible upper bound 
and secondly, it relies on the observed spectrum being white. Unfortunately, even 
for millisecond pulsars, it is unusual for the timing residuals not to be affected by 
unexplained phenomena, such as pulsar "timing noise" (see e.g. Hobbs, Lyne & Kramer 
2006). A new method has now been developed and will be described in a subsequent 
paper. Our new method is applicable to any pulsar timing data-set and will lead to 
significantly improved bounds on the GW background amplitude. 

3. Detecting a background 

Jenet et al. (2005) developed a technique for making a definitive detection of an 
isotropic, stochastic, background of GWs by searching for correlated signals between 
pulsar data-sets. For each pair of pulsars, the zero-lag correlation between the respective 
timing residuals can be calculated. The expected correlation as a function of angle 
between the pulsars for timing residuals dominated by a GW background was first 
calculated by Hellings & Downs (1983) and is shown in Figure |lfj|. Jenet et al. 
(2005) showed, for expected amplitudes of a GW background, that this signal could 
be unambiguously detected if 20 or more pulsars were observed over a period of 5 years 
each with an rms timing residual of 100 — 500 ns. 

We have recently generalised the initial work of Jenet et al. (2005) to produce 
routines that can be applied to the observed pulsar timing residuals (with the actual 
sampling, data lengths etc.). This new technique also increases the significance of 
any possible detection, by optimally pre-whitening the data-sets. An example of the 
improvement possible by pre-whitening the timing residuals is shown in Figure [5] where 
the significance of detection is plotted versus the GW background amplitude for the 

|| This correlation curve assumes the general theory of relativity. We are calculating the equivalent 
curves for other theories of gravity. These curves and their implications will be discussed elsewhere. 
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Figure 4. The expected correlation in the timing residuals of pairs of pulsars as a 
function of angular separation for an isotropic GW background. 

PPTA goal of timing 20 pulsars with 100 ns timing precision over a period of five 
years. Without any pre-filtering of the data we will only achieve a maximum detection 
significance of S ~ 3.5 (corresponding to the sensitivity labelled as the "PPTA limit" 
in Figure [T]). Using optimal whitening schemes this sensitivity can be significantly 
enhanced. Initial results suggest that if we can reduce our current timing precision for 
all pulsars by a factor of ~ 2 we will be sensitive, within five years, to the maximum 
amplitude predicted by models of an isotropic stochastic GW background caused by 
coalescing black-holes in the centres of galaxies^- 

4. Practical issues 

Pulsar timing is affected by the stability of terrestrial clocks, ephemeris errors and the 
pulsars themselves. The goal of the PPTA project, described above, to obtain data-sets 
for 20 pulsars with approximately two-weekly sampling with rms timing residuals of 
a ~ 100 ns over five years, is challenging. To date we have a few pulsars where we can 
obtain such precision. For instance, van Straten et al. (2001) obtained a = 130 ns over 
40 months of observing PSR J0437— 4715. For the majority of our pulsars o ~ 1/is. We 
have many possibilities for decreasing these timing residuals, including: 

• correcting the data-sets for dispersion measure variations: You et al. (2007) showed 
that for PSR J1939+2134, a decreased from 290 ns to 190 ns after correction. 

• improving instrumentation: We have recently developed new digital filterbank 
systems for the Parkes telescope. The recently commissioned system provides 1 GHz 

% Note that this limit is similar to the current limits described in §2.3 that were obtained using a few 
long data-sets. We emphasise that, even though limits can be placed using a few data-sets, detecting 
a GW background requires observations of ~20 pulsars. 
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Figure 5. The sensitivity of the completed PPTA project to a GW background of 
specified amplitude. The three curves represent no pre-filtering of the data (lowest), 
using an optimal low-pass filtering technique (middle) and optimally pre-whitening the 
data-sets (highest). 

of bandwidth at an observing frequency of 3 GHz. We are also developing a new 
wide-bandwidth coherent dedispersion system. 

• new processing techniques: We are exploring new methods to calibrate our data-sets 
which should lead to at least a factor of 2 improvement in our rms residuals. We 
have recently produced a new pulsar timing package, tempo2 (Hobbs, Edwards & 
Manchester 2006), that is accurate (for known physics) at the Ins level. 

5. Future possibilities 

The sensitivity of a pulsar timing array to a stochastic background is proportional to 
the average timing precision, the square root of the number of observations and to 
the number of pulsars in the array. This paper has concentrated on the Parkes pulsar 
timing array project. However, significant improvements will be made to the timing 
array sensitivity to GW sources by combining our data-sets with Northern Hemisphere 
pulsar observations. The North American pulsar timing array project (NanoGrav) has 
been observing a sample of pulsars using the Arecibo and GreenBank observatories. 
These observations provide more than 20 years of timing PSRs B1855+09 and B1937+21 
and were used to provide the first stringent limits on the stochastic gravitational wave 
background (e.g. Kaspi et al. 1994) and to rule out the postulated binary black hole 
system in the radio galaxy 3C66B (Jenet et al. 2004). The European Pulsar Timing 
Array (EPTA) team have recently started to obtain data on a large sample of pulsars 
using the four major European telescopes at Jodrell Bank, Effelsberg, Nancay and 
Westerbork. We are currently developing techniques to combine observations from 
different observatories and will soon produce new limits on the existence of a GW 
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background using all available data. 

Figured] shows the sensitivity that should be achievable using the Square Kilometer 
Array (SKA) telescope hoped to be built by the year 2020. This sensitivity curve 
assumes an SKA timing array project that could observe 100 millisecond pulsars with 
rms timing residuals around 50 ns for a period of 10 yr. Provided intrinsic pulsar timing 
noise does not dominate the effect of the stochastic GW background then this limit 
corresponds to a detection limit at 3nHz of Q gw ~ 10~ 13 . Pre-whitening procedures 
have the potential to greatly decrease this limit. 

6. Conclusion 

It is possible that gravitational waves will be detected within the next decade by world- 
wide pulsar timing array projects. As a by-product of these investigations stringent 
checks will also be placed on terrestrial time standards and the solar system ephemeris. 
Regular dual-frequency observations of multiple pulsars will also provide valuable 
information about the interstellar medium. Using a pulsar array as a gravitational 
wave detector is complimentary to other searches that are attempting to detect much 
higher-frequency gravitational waves. 
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